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I study the H C 2 transition within the Ginzburg-Landau model, with m-component order pa- 
rameter tpi. I find a renormalized fixed point free energy, exact in m — > oo limit, suggestive of a 
2nd-order transition in contrast to a general belief of a lst-order transition. The thermal fluctuations 
for H ^ force one to consider an infinite set of marginally relevant operators for d < d uc = 6.1 find 
di c = 4, predicting that the ODLRO does not survive thermal fluctuations in d = 2, 3. The result 
is a solution to a critical fixed point that was found to be inaccessible within e = 6 — d-expansion, 
previously considered in E.Brezin, D.R.Nelson, A.Thiaville, Phys.Rev.B 31, 7124 (1985), and was 
interpreted as a lst-order transition. 
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The nature of the Normal-to-Superconducting (NS) 
transition is substantially more challenging than the neu- 
tral superfluid transition and has remained controversial 
for many years. While there is now finally a general con- 
sensus on the nature of the zero field NS transition [1-3], 
the finite field NS transition is much less understood. As 
discovered by Abrikosov [4], in mean-field theory, upon 
lowering T the NS transition takes place at H C 2 (T C (H)) 
and the superconducting order parameter develops a tri- 
angular array of flux- line vortices, a transition believed 
to be describable by the Ginzburg-Landau free energy 
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where r ~ (T — T c )/T c , q = 2e/hc, and [i is magnetic 
permeability of the normal metal and is close to unity. 
Although it is generally believed that in mean-field the 
transition is of second order, even this theory is not on a 
completely firm ground. The issue is that at the transi- 
tion two apparently unrelated symmetries are supposed 
to be simultaneously broken: the translational symmetry 
is broken by the formation of the Abrikosov lattice ac- 
companied by breaking of U(l) gauge symmetry leading 
to phase coherence and superconductivity. While it is 
easy to see the diverging uniform susceptibility (k = 0), 
associated with the formation of phase coherence, no cor- 
responding divergent susceptibility for the formation of a 
lattice (much less a specific triangular lattice) had been 
demonstrated. Following Abrikosov, what is generally 
demonstrated is that a periodic solution (e.g. triangu- 
lar) is lower in energy, globally, than a lattice of another 
symmetry (e.g. the square one). 

The gaussian fluctuation corrections only to the spe- 
cific heat around H C 2 (valid only outside critical re- 
gion) were first studied in the pioneering work by Lee 
and Shenoy [5], followed by a more sophisticated self- 



consistent method and expansions in the range of the 
nonlinear interaction [6]. These results for the specific 
heat are consistent with my more general and controlled 
analysis of the full critical behavior. 

Although a form of freezing into a 3d lattice takes 
place, as was first emphasized by Brezin, Nelson and Thi- 
aville (BNT) [7], Landau's general argument for the lst- 
order freezing transition is circumvented by the fact that 
in mean-field both of the symmetries are broken simulta- 
neously and the amplitude of ip is small, preventing the 
density from jumping discontinuously. Since it is likely 
that this will no longer be the case, once thermal fluc- 
tuations are taken into account, (the formation of the 
lattice will be depressed below H 7 ^ 1 , where amplitude 
of ip forms), it is generally believed that the transition 
will become lst-order. [7] 

Thermal fluctuations were first studied by BNT [7] in 
e = 6 — d expansion about the upper critical dimen- 
sion d uc = 6. They found that fluctuations force one 
to consider an infinite set of marginally relevant quar- 
tic operators and derived a one-loop functional renor- 
malization group recursion equation for the quartic cou- 
pling function. However, numerically analyzing these 
integro-differential equation, they were not able to find a 
fixed point function within the e-expansion. Consistent 
with their physical argument, BNT interpreted the run- 
away rg flows toward a negative coupling function as a 
fluctuation-driven lst-order transition. 

In clean samples the experiments observe a lst-order 
melting transition of the Abrikosov vortex lattice [8] , con- 
sistent with the conclusions of BNT. However, more re- 
cent experiments have found that this line of the lst- 
order transition terminates at a critical end point, at a 
well-defined magnetic field, giving way to a 2nd-order 
transition [9]. 

In this Letter, I study the H C 2 NS transition by gen- 
eralizing the GL theory to a large number m of complex 
components of the superconducting order parameter %pi 
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and solve it in the m — > oo limit in arbitrary dimension- 
ality. I find a 2nd-order transition in contrast to a gen- 
eral belief of a lst-order transition, which, however only 
survives for d > di c = 4 and therefore predict that the su- 
perconducting phase coherence does not survive thermal 
fluctuations in physical dimensionality (d = 2,3). 

Near H c2 , considering strongly type-II superconduc- 
tors allows me to ignore gauge field fluctuations. I 
take A = ^H(—y,x;0±), as the symmetric gauge d- 
dimensional generalization of H field perpendicular to 
the x-y plane. I solve the problem in the lowest Lan- 
dau level (LLL) approximation valid near H c2 . In this 
case the quadratic part of the Hamiltonian Eq.l can be 
diagonalized by ipi(z, z; fj_) = 4>i(z;r)e~\ z \ Z 4 ' , where 
4>i(z;r) is an arbitrary m-component analytic function 
of z = x + iy (not to be confused with the z coordinate), 
f± denotes the d — 2-dimensional space perpendicular to 
the x-y plane defined by the H field, and I = 1/y/qH is 
the magnetic length. Substituting this form for ipi into 
Eq.l I obtain, 

F[<t>i] = J d d -r ± d 2 z ie -^^ 212 (|Vx&| 2 +t |&| 2 ) + 
\ j d 2 z 2 e~^ 2 / 212 g (\ Zl - z 2 \)\Uzi,f±)\ 2 \U^,r ± )\ 



expect a finite T c NS transition. 

I now follow a standard large m treatment [10], which 
surprisingly, can be easily generalized to the problem at 
hand. Via a Hubbard-Stratanovich transformation I in- 
troduce an auxiliary field x( z >r±) that linearly couples to 
\(t>i(zi,r±)\ 2 and has a bare propagator g (\zi — z 2 \) with 
the inverse defined by J d 2 zg^(\zi — z\)g a (\z — z 2 \) = 
(z\ — z 2 ). Anticipating a breaking of symmetry along 
one of the i-directions, I take 4>i = and inte- 

grate over the m — 1 complex (Goldstone modes to-be) 
components of 7r Q . For g {\z\ — z 2 \) of order 1/m, in the 
large m limit, the fluctuations in the fields a and x can 
be ignored (they lead to 1/m corrections) and I obtain 
the exact fluctuation corrected effective free energy 
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where t Q = r Q (T) + qH oc (H - H c2 (T)) ex. (T — T C (H)) 
changes sign at the mean-field theory (MFT) transition 
and the interaction coupling constant g has been gen- 
eralized from g 8^ 2 \z\ — z 2 ) to an arbitrary short-range 
function g {\z\ — z 2 \) evolution of which to a fixed point 
function (due to thermal renormalization) is the main re- 
sult of this work. The fact that one must keep track of 
the whole function <?(|;z|), rather than expanding can be 
seen as follows. In order to reach the NS transition fixed 
point, while r± is treated as usual distances rescaling un- 
der dilation by b as, r± — > r±b, z must remain unrescaled 
z — > z, in order to preserve under rescaling the form of 
the quadratic part of the free energy. If z were rescaled 
with any power of b, the e~l z l I 21 would reduce the free 
energy to zero, as b — > oo, clearly an unsatisfactory sit- 
uation. The nondimensionality of z therefore leads to 
an infinite set of quartic coupling constants, labeled by z 
and encoded into the g(\z\) coupling function. Physically 
this behavior is due to the fact that while the correlation 
length £j_ diverges (I rescale r± to keep up with this 
divergence), the q = wavevector correlations in the x- 
y plane grow but are eventually cutoff by the magnetic 
length I. All the long- wavelength modes labeled by z are 
equally important in contributing to the divergences in 
k± — > susceptibility and other correlation functions, 
and therefore must be treated on equal footing. Rescal- 
ing of z — > zb would amount to looking at the theory 
characterized by a vanishing magnetic length I — > l/b or 
equivalently divergent H field, a regime in which I do not 
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To make analytical progress in evaluating the free energy, 
I will make the following choice for the low energy ansatz, 

a(z,r±) = a G a G (z) = a G JJ(z - z t ) , 

i 

x(z,f±)=Xo, (4) 

and now explain the physically motivation behind this 
choice. Obviously taking a(z,f±) to be uniform along 
fx lowers the free energy. Since I want to describe the 
transition into an Abrikosov-type of phase, i.e. a phase of 
periodic arrangement of vortices characterized by a G (z) 
(probably a triangular lattice given by a global minimiza- 
tion procedure), I take a(z) to be proportional to this a 
priori assumed low T phase. The complex amplitude co- 
efficient a G characterizes the level of order in this phase. 
Since the square of the superconducting order parame- 
ter |-0(2;)| 2 , is essentially constant (aside from a lattice of 
zeros), J d 2 z\a(z)\ 2 e~\ z \ I 21 = Aaa 2 G = Ap s is extensive 
with the x-y area A (a is the reduction factor in the su- 
perfluid density p° s due to zeros). Since x( z ^±) couples 
to \<t(z, r±)\ 2 , with equivalent thermal averages, it is a 
good approximation to take it to be a constant, as I have 
done in Eq.4. Certainly the nonuniformity coming from 
the zeros in Xo (z) will lead to numerical corrections but it 
is unlikely that they will result in any modification of the 
universal and qualitative properties of the NS transition. 

Another important point is that I am obviously as- 
suming the form of the low-T phase, drawing on the 
knowledge of a lattice solution. This is in the spirit of 
Landau's theory of freezing where the effective theory of 
the transition is for the slowly varying coefficient p G of 
the particular Fourier component of the full density at 
which the ordering is to occur. As in the theory of freez- 
ing, it is an exceptionally difficult task to analytically 
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demonstrate a local instability to this finite wavevector 
G ordering, which is brought about by the repulsive in- 
teractions. For the theory of freezing one does a separate 
difficult calculation within the disordered (liquid) phase 
to show that interactions lead to a peak in the struc- 
ture function, approximately at the inverse of the inter- 
particle separation, this peak being the precursor of the 
inevitable Bragg peak. The analogous task of showing 
(within the GL theory) a susceptibility to order at a fi- 
nite in-plane wavevector, into for example the Abrikosov 
lattice, has not, to my knowledge, been performed an- 
alytically, even at the MFT level. Instead, one simply 
assumes a periodic solution and picks the global lowest 
energy one. Since the fluctuation-corrected free energy 
in Eq.3 is significantly different from the MFT one, it is 
possible, that upon minimization, it will lead to a lattice 
of a different form from the triangular one. 

Having "apologized" for my choice of the ansatz, I now 
minimize F in Eq.3 with respect to cto and Xo, obtaining 
the saddle point equations 



H field quenching of the kinetic energy in the x-y plane, 
perpendicular to it. 

I now turn to the main focus of my work, that is, the 
effective free energy that describes the fluctuations about 
this large m solution. In the large m limit, the renormal- 
ized quartic interaction gn(\zi — Zi\) is the quantity that 
characterizes the fixed point of the NS transition. One 
way of computing it is to go back to the free energy in 
Eq.2 and to simply pertubatively resum all the loop dia- 
grams to lowest order in 1/m, which renormalize the bare 
quartic interaction g (\zi — Zi\). Equivalently, given the 
formalism that I have set up above, I can simply com- 
pute the renormalization of the quadratic term in \ by 
expanding the free energy in Eq.3 to quadratic order in 
fluctuations of x about the saddle point value Xo- Either 
approach gives, 
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where q is a wavevector in the x-y plane, II(g, k±, t) is the 
Fourier-transform of the polarization bubble with respect 
to complex z (x-y) coordinate and 



where b = 2/(ma) J d 2 zg 1 (\z\), and I defined the cor- 
relation length £ by £~ 2 = t Q + Xo- For T < T C (H), 
aa > 0, and Eq.5 implies a divergent correlation length 
(£ -2 = 0) in the Abrikosov phase. Eq.6 then leads to the 
order parameter aa, vanishing as a power-law with the 
reduced temperature, as the T C (H) is approached from 
below 
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to thermal fluctuations. I note that for d = d± + 2 < 4 
the thermal suppression of T c is divergent and therefore 
identify the lower-critical dimension as d; c = 4. This is 
consistent with conclusions about d/ c by Moore based on 
a completely different arguments applied in the ordered 
Abrikosov phase [11]. Eq.7 implies a vanishing of the su- 
perfluid density as T C (H) is approached from below as 
p s ~ | QG |2 _ (T C (H) -T). For T > T C (H) the solution 
is characterized by a finite £. Eq.5 then implies aa = 
and the state is normal. Eq.6 then describes how the cor- 
relation length £ diverges as T C (H) is approached from 
above 

£ ~ (T — T C (H))~" , (8) 

where for d > d uc = 6 the MFT is accurate and v = \, 
while for di c < d < d uc thermal fluctuations are diver- 
gent and lead to v = l/(d — 4). These are, as expected, 
the usual large m exponents in the d — 2 dimensional the- 
ory, with the dimensional reduction occuring due to the 



where in Eq.12 I(k±,t) is evaluated at T c (t = 0), and 
c(d) = r(2 - d ± /2)T(dj2 - l) 2 /T(d ± - 2)/(47r)^/ 2 , 
diverging as 1/(6 — d) near d uc , consistent with the e- 
expansion of BNT. The interaction function gn(q, k±,t = 
0), at T c , is displayed in the inset of Fig.l as a function 
of q, for a fixed small k±, large m, and for a choice of 
is the shift in T C (H) due short-range microscopic GL interaction g (q) = e~ a q / 2 , 
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with a 

q = G « Mmil^y/lhik 



1. It shows a peak at a finite wavevector 
a -1 }, as one might expect 



for a continuous freezing transition, with a peak possi- 
bly being a precursor of the eventual Bragg singulari- 
ties, not captured by the large m theory. In real space 
gR( r xy, k±, t = 0) exhibits oscillations with a 1/r 2 ^ power 
law fall off, reflecting correlations building up in the xy 
plane. Eqs.9-12 contain the information about the MFT 
fixed point (valid for d > 6), the new Heisenberg-type 
fixed point, characterizing the NS transition in the pres- 
ence of thermal fluctuations, and about the crossover be- 
tween them as a function of the reduced temperature t 
and length scale k±, labeled by the x-y wavevector q. 
Right at T C (H), for fixed q, U(q,k±,0) ~ fc x € diverges 
at long scales and MFT crosses over (with a crossover 
exponent <j> = — e) to the new fixed point summarized 
by Eqs.5-8, and with the effective interaction for fluctu- 
ations given by a completely universal function, (actu- 
ally a distribution) , independent of the bare interaction 
g (\zi - 22 |) 
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where g = 4tt1 2 /m/c(d). This form of gn is valid for ar- 
bitrary d, but for 6 — d = e « 1 gives the fixed point 
function of BNT rg equations, in to — > oo limit [12]. 



>g(r xy ,k±,0) 



g(q,k±,0) 




FIG. 1. gR{r xy ,k±,t = 0) at T c for fixed k± exhibits den- 
sity correlations in the xy plane manifested in oscillations at 
a scale G _1 ; the inset shows gn(q,k_i_,t = 0). 

The fluctuations near H C 2 at the NS critical point are 
described by the effective free energy of same form as 
Eq.2, but with g a replaced by the universal function gn 



from Eq.14. Because x couples to 



it is easy to show 



that its renormalized propagator gn is also the superfluid 
density susceptibility. This allows me to extract the spe- 
cific heat a = — (6 — d)/(d — 4) and the correlation length 
v = l/(d — 4) exponents, consistent with the value of 
v obtained from the saddle point equations. Note that 
because of the unusual scaling of lengths at this critical 
point, discussed in the introduction, the standard rela- 
tion a = 2 — dv is replaced by an easily derivable relation 
appropriate for this problem, a = 2 — (d — 2)v. The 
screening by fluctuations has removed the MFT specific 
heat divergence (a m f t = (6 — d)/2 > 0) and replaced it 
by a nondivergent nonanalyticity in k± with a < 0. This 
suggests by a variant of the Harris criterion that the in- 
finitesimal short-range point disorder is irrelevant for this 
finite field NS transition [12]. Similar arguments suggest 
that the transition is destabilized by the long-range dis- 
order such as twin planes and the artificially introduced 
columnar defects [13]. 

Two of the important remaining problems are the 
demonstration of a diverging susceptibility to order into 
a lattice of vortices (i.e. divergent £ xy at q ~ G), and the 
computation of the 1/m-corrections. Preliminary calcu- 
lations indicate that these problems might be related, 
and that the renormalized interaction in Eq.14 can lead 
to an instability to a vortex lattice, when the 1/to cor- 



rections to the effective free energy are computed [13]. 

In summary, I have presented a large to theory of a 
2nd-order H C 2 NS transition for arbitrary d and found 
the fixed point effective free energy describing the tran- 
sition. Since the physical dimensionalities d = 2,3 are 
below the di c = 4, the transition does not survive ther- 
mal fluctuations in the limit of large to. It is, however, 
possible that for the physical value of m = 1, that the 
true lower-critical dimension will be reduced below d = 3, 
allowing the finite-T, continuous NS transition to persist. 

Note added: After this work was submitted for publi- 
cation I became aware of a interesting but unfortunately 
relatively unknown paper by I. Affleck and E. Brezin, 
Nucl. Phys. B 257, 451 (1985) that also treats H c2 tran- 
sition iiiHUoo limit. In contrast to my work, where 
I focus on the fluctuation part of the renormalized free 
energy, that work studies only the m = oo saddle point 
equations to the constant part of the free energy (my 
Eq.3). Affleck, et al. argue that there is no solution to 
the saddle point equations near H C 2, interpreting it as a 
lst-order transition, in contrast to my conclusion. We 
have not been able to reconcile a possible absence of a 
saddle point solution and my unambiguous finding of the 
fixed point free energy (determined by Eq.14), usually 
indicative of a 2nd order transition. 

I thank M. Moore for many enlightening discussions. I 
acknowledge the support by the NSF (DMR 94-16926), 
through the Science and Technology Center for Super- 
conductivity. I am grateful to The Institute for The- 
oretical Physics at UCSB, where this work was done, 
for their hospitality, and their support under NSF Grant 
No. PHY89-04035. 



[1] B. I. Halperin, T. C. Lubensky and S. K. Ma, 

Phys. Rev. Lett. 32, 292 (1974). 
[2] C. Dasgupta and B. I. Halperin, Phys. Rev. Lett. 47, 

1556 (1981) 

[3] L. Radzihovsky, Euro. Phys. Lett. 29, 227 (1995). 

[4] A. A. Abrikosov, Zh. Eskp. Teor. Fiz 32, 1442 (1957) 

[Sov. Phys. - JETP 5, 1174 (1957)]. 
[5] P. Lee and S. Shenoy, Phys. Rev. Lett. 28, 1025 (1972). 
[6] A. Bray, Phys. Rev. B 9, 4752 (1974); S. Hikami and A. 

Fujita, Phys. Rev. B 41, 6379 (1990). 
[7] E. Brezin, D. R. Nelson and A. Thiaville, Phys. Rev. B 

31, 7124 (1985). 
[8] H. Safar et al, Phys. Rev. Lett. 69, 824 (1992). 
[9] H. Safar et al, Phys. Rev. Lett. 70, 3800 (1993). 
[10] J. Zinn- Justin Quantum Field Theory and Critical Phe- 
nomena (Oxford University Press, NY, 1989). 
[11] M. A. Moore, Phys. Rev. B 45, 7336 (1992). 
[12] M. A. Moore, T. Newman and L. Radzihovsky, unpub- 
lished. 

[13] L. Radzihovsky, unpublished. 



4 



